Permanent magnetic lattices for ultracold atoms 
and quantum degenerate gases 

Saeed Ghanbari, Tien D Kieu, Andrei Sidorov and Peter 
Hannaford 

Centre for Atom Optics and Ultrafast Spectroscopy and 
ARC Centre of Excellence for Quantum Atom Optics 
Swinburne University of Technology, Melbourne, Australia 3122 

Abstract. We propose the use of periodic arrays of permanent magnetic films 
for producing magnetic lattices of microtraps for confining, manipulating and 
controlling small clouds of ultracold atoms and quantum degenerate gases. Using 
analytical expressions and numerical calculations we show that periodic arrays 
of magnetic films can produce one-dimensional (ID) and two-dimensional (2D) 
magnetic lattices with non-zero potential minima, allowing ultracold atoms to be 
trapped without losses due to spin flips. In particular, we show that two crossed 
layers of periodic arrays of parallel rectangular magnets plus bias fields, or a single 
layer of periodic arrays of square-shaped magnets with three different thicknesses 
plus bias fields, can produce 2D magnetic lattices of microtraps having non- 
zero potential minima and controllable trap depth. For arrays with micron-scale 
periodicity, the magnetic microtraps can have very large trap depths (~0.5 mK 
for the realistic parameters chosen for the 2D lattice) and very tight confinement. 



1. Introduction 

Periodic optical lattices produced by the interference of intersecting laser beams have 
been used extensively in recent years to confine, manipulate and control small clouds 
of ultracold atoms and Bose-Einstein condensates pQ. Such lattices are ideal tools 
for performing fundamental quantum physics experiments such as studies of low- 
dimensional quantum gases [2] and quantum tunnelling experiments including the 
BEC superfiuid to Mott insulator quantum phase transition [3]. Optical lattices also 
have potential application in quantum information processing since they may provide 
storage registers for qubits based on neutral atoms |S] . 

An alternative approach for producing periodic lattices for ultracold atoms is to 
use the magnetic potentials of periodic arrays of magnetic microtraps. Simple, one- 
dimensional (ID) magnetic lattices consisting of arrays of 2D traps or waveguides 
have been proposed [H] and constructed using current-carrying wires [Jj or permanent 
magnetic structures |B1 El Ell on 'atom chips', and two-dimensional (2D) lattices of 
magnetic microtraps produced by crossed arrays of current-carrying wires have been 
proposed [TTIIT!?] . 

We have recently developed technology for producing high-quality magnetic 
microstructures based on permanent, perpendicularly magnetised magneto-optical 
Tb e Gd w Feg,oCo4 films These magnetic microstructures have been used to 

construct periodic grooved magnetic mirrors for ultracold atoms |13l 114] , which in 
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the presence of bias magnetic fields can be transformed into a ID magnetic lattice 
of 2D traps, and to produce atom chips for ultracold atoms and Bose-Einstein 
condensates ^1- The Tb 6 Gdi FesoCo4 magneto-optical films exhibit excellent 
magnetic properties for atom optics applications: they can be deposited with high 
perpendicular magnetic anisotropy; they have excellent magnetic homogeneity; and 
they can have high permanent magnetisation (AirMz ~ 3.8 kG), large coercivity 
(H c ^3 kOe) and relatively high Curie temperature (~300°C). 

In this paper we propose the use of periodic arrays of permanent magnetic films for 
producing magnetic lattices of microtraps for confining, manipulating and controlling 
small clouds of ultracold atoms and quantum degenerate gases. Using analytical 
expressions and numerical calculations we show that it is possible to produce ID and 
2D permanent magnetic lattices with non-zero potential minima in which ultracold 
atoms prepared in low magnetic field-seeking states can be trapped without losses due 
to Majorana spin flips. In particular, we show that two crossed separated layers of 
periodic arrays of parallel rectangular magnets plus bias magnetic fields, or a single 
layer of periodic arrays of square-shaped magnets having three different thicknesses 
plus bias fields, can produce 2D magnetic lattices of microtraps with non-zero potential 
minima and large and controllable trap depth and high trap frequency. 

Magnetic lattices based on permanent magnetic films have potential advantages 
over optical lattices or magnetic lattices based on current-carrying wires that make 
them attractive for atom optics applications and compact integrated devices. They 
do not involve (high intensity) laser beams or any beam alignment, and there is 
no light scattering or decoherence due to spontaneous emission. They can produce 
highly stable reproducible potential wells with low technical noise. They can be 
produced with large trap depth and large magnetic field curvature, leading to 
very high trap frequencies, without heat dissipation, in contrast to current-carrying 
microwire devices. Using modern microtechnology, permanent magnetic lattices may 
be fabricated with a wide range of periods from about 100 fim down to about 1 \xm 
and they can in principle involve variable spacing between the lattice sites or complex 
potential shapes at each lattice site. Finally, in magnetic lattices, only atoms in 
low magnetic field-seeking states are trapped and it should be possible to perform 
radiofrequency evaporative cooling in situ in the lattice, thereby allowing the study 
of very low temperature phenomena in lattices. 

2. Analytical expressions for infinite magnetic lattices 

2.1. Single infinite periodic array of magnets with bias fields 

We consider first the simple case of a single infinite periodic array of parallel, 
rectangular, long magnets of thickness t, with periodicity a along the y-direction, 
perpendicular magnetisation M z , and uniform bias magnetic fields B\ x , B\ y and B\ z 
along the x-, y- and z-directions [figure IJa) -(c)]. 

The components of the magnetic field at distance z from the bottom surface of 
the array of magnets [figure ^b)] can be written as (using the results given in 
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Figure 1. (a-c) Single periodic array of parallel, rectangular magnets with 
perpendicular magnetization. (d-f) Contour plots of the magnitude of the 
magnetic field in the central region in the yoz plane without bias fields (d), 
with a bias field Bi y = —15 G along the y-direction (e), and with bias fields 
B\ x = —20 G, Bi y = — 15 G and B\ z = —0.09 G along the x-, y- and ^-directions. 
For this calculation, the number of magnets n r = 1001, t = 0.05 fim, a = 1 fim, 
l x = 1000.5 fim, and 4nM z = 3.8 kG. The spacing between the contour lines is 
7 G. 



B z = B [(l - c - kt )e- k ^-^ s+t ^ cos(fcy) 

- i(l - e ~3kt^-3k[ z -{s+t)] cos ( 3fcy ) + ...]+ Blz (lc) 
o 

where the decay constant k = 2ir/a, B = 4M Z (Gaussian units), and the factors 
(1 — e~ fet ), (1 — e~ 3kt ), ■ ■ ■ account for the finite thickness t of the magnets. 

For distances from the surface which are large compared with a /An, the higher 
order spatial harmonics in Ijl&t and Ijlcfl will have decayed away, and (|la|) - ((Tc|) reduce 
to 

B x = B lx (2 a) 

By = B 0y sm(ky)e- kz + B ly (26) 
B z = B 0y cos(ky)e~ kz + B lz (2c) 
where B 0y = B (l - e~ kt )e kt . The magnitude of the magnetic field is then given by 
B(y,z) = [Bl x + Bl y + B 2 lz 

+ 2[B 0y B ly sin(fcy) + B 0y B u cos(fcy)] e - fc2 + Bl y e~ 2kz } 2 (3) 

The effect of the bias fields B\ y and B\ z is essentially the same; so we set B\ z = 0. 
This results in a ID periodic magnetic lattice of 2D magnetic traps with non-zero 
potential minima, given by 

Bmin = | -Biz | (4) 
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Umin = ( + 7 j " (5a) 

(56) 

where n y = 0, ±1, ±2, • • • represents the trap number in the y-direction. The barrier 
heights in the y- and z-directions are given by 

AB* = Bl ax - B mm = (Bf x + 4Bf y )i - \B lx \ (6a) 
AB Z = B* max - B mm = [B\ x + B\ y )^ - \B lx \ (66) 

The curvatures of the magnetic field at the centre of the 2D magnetic traps and the 
trap frequencies in the y- and z-directions are given by 



d 2 B _ d 2 B _ 4tt 2 B\ y 



Qy2 Q z 2 a 2 \Bi x 

2tt / m F g F fi E 



a V rn\Bix\ 



(7) 



\Bly\ (8) 



where m F is the magnetic quantum number of the hyperfine state F, g F is the Lande 
g-factor, [i B is the Bohr magneton and m is the atomic mass. Equations (|5^l) and (JSJl 
are the same as given in [S]. 

2.1.1. No bias fields For the case of no bias fields, © reduces to 

B(z) - B 0y e~ kz (9) 

Under these conditions the magnitude of the magnetic field falls off exponentially with 
distance z from the surface. This configuration represents the familiar case of a grooved 
magnetic mirror for slowly moving atoms in low magnetic field-seeking states |14| . 

2.1.2. Single bias field B\ y For the case of a single bias field Bi y , Q becomes 

B(y,z) = [B\ y + 2B 0y B ly sm(ky)e~ kz + ^e^}* (10) 

The magnitude of the magnetic field develops corrugations with period a in the y- 
direction, and 2D magnetic traps with B m i n — appear in the potentials above the 
surface. This configuration, involving a single bias field, may be used as a spatial 
diffraction grating for slowly moving atoms ^1 El El 1201 > but is not useful as a 
magnetic lattice for ultracold atoms because of the zero potential minima. 
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Figure 2. Periodic array consisting of two crossed layers of parallel, rectangular 
permanent magnets with perpendicular magnetisation. In (d) the locations of the 
central minima in the xoy plane are shown for a symmetrical magnetic lattice 
with bias fields B lx = -4.08 G, B ly = -6.05 G and B lz = -0.69 G and other 
parameters as given in table 1, column 3. 



Table 1. Input parameters for (1) two crossed layers of periodic arrays of 
rectangular magnets (figure |5J, and (2) a single layer of square magnets with 
three thicknesses (figure 0. 



Parameter 


Definition 


Configuration (1) 


Configuration (2) 


n r or n sq 


Number of rectangular 


n r = 1001 


n 3q = 401 




magnets or square magnets 








in x- or {/-direction 






a (fim) 


Period of magnetic lattice 


1.000 


1.000 


lx = ly (pm) 


Length of magnets or 


1000.5 


200.5 




magnetic array along x or y 






s (p.m) 


Separation of magnetic layers 


0.100 




ti (/Ltm) 


Thickness of magnetic 


0.322 


0.120 




film (first) 






ti (nm) 


Thickness of magnetic 


0.083 


0.100 




film (second) 






t 3 (/Ltm) 


Thickness of magnetic 




0.220 




film (third) 






4ttM z (G) 


Magnetisation along z 


3800 


3800 


B lx (G) 


Bias magnetic field along x 


-4.08 


-5.00 


Bly (G) 


Bias magnetic field along y 


-6.05 


-4.22 


B u (G) 


Bias magnetic field along z 


-0.69 


-1.87 
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2. 2. Two crossed layers of infinite periodic arrays of magnets with bias fields 

We now consider a configuration consisting of two crossed, separated, infinite periodic 
arrays of parallel, rectangular, long magnets with perpendicular magnetization M z and 
uniform bias fields B\ Xl B\ y and B\ z along the x-, y- and z-directions (figure |3J). The 
bottom array has periodicity a along the y-direction and the top array has periodicity 
a along the x-direction. The two arrays are separated by a distance s and the magnets 
in the bottom array have thickness t\ while those in the top array have thickness ti 
[see figure |2Jb)]. This configuration of crossed permanent magnets with bias fields 
has similarities to the crossed current-carrying wire configurations proposed by Yin et 
al and by Grabowski and Pfau • The components of the magnetic field for 
distances from the top surface which are large compared with a/4ir are given by 

B x = B Qx sm{kx)e~ kz + B lx (11a) 

B y = B Qy sm(ky)e- kz + B ly {lib) 

B z = [B 0x cos(fcx) + B 0y cos(ky)]e- kz + B u (11c) 



where 



B 0x = B Q (1 - e - kt2 )e kis+tl+t2) (12a) 
B 0y -B (l-e- fctl )e fetl (126) 
The magnitude of the magnetic field above the crossed magnetic arrays is then 
B(x,y,z) = {B 2 lx + Bf y + Bl 

+2[B 0x Bi x sin(kx) + Bo y B\ y sin(fcy) 
+B Qx B lz cos{kx) + B 0y B lz cos(ky)]e- kz 

+ [Bq x + B 2 Q y + 2B Qx B 0y cos(kx) cos(ky)]e~ 2kz } 2 (13) 



2.2.1. Two bias fields B\ x and B\ y For the case of two bias fields B\ x and B\ y , (| 1 31) 
becomes 



?2 



B(x,y,z) = {B 2 lx + Bl 

+2[B 0x B lx sm(kx) + B 0y B ly sm(ky)}e~ kz 

+ [B 2 Qx + Bl y + 2B 0x B 0y cos(kx) cos(fcy)]e- 2 ' £Z } 2 (14) 

This configuration results in a 2D periodic lattice of magnetic traps with non-zero 
potential minima given by 

\Br\^Bi — Br\',,B-\ ^ 



or 
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2.2.2. Symmetrical magnetic lattice with bias fields B lx and B Xy To produce a 
symmetrical magnetic lattice the amplitude of the oscillating magnetic field produced 
by the bottom magnetic array in the y-direction needs to equal the amplitude of the 
oscillating magnetic field produced by the top array in the ^-direction. To satisfy this 
condition we impose the constraint 

Bq x Bix = B 0y Bi y (16a) 
B ly = c B lx (16b) 

where 

is a dimensionless constant which only involves geometrical constants a, s, t\ and i 2 of 
the magnetic array. The magnetic traps then have non-zero potential minima given 
by 

Brnin = C\\Bi x \ (18) 

which are located at 

Xmin = [n x + ^ a, n x = 0,±1,±2, ••• (19a) 

J)a, % = 0,±1,±2,--- (196) 
a ( c 2 B 0x \ 

Zmin = — m ( T^j-T J ( 19C ) 



2tt V LB 



where c\ and c 2 are dimensionless constants which may be expressed in terms of the 
constant c = B 0x /B 0y 

11 — r 2 \ 

d = 11 C °\ (20a) 

c 2 4(l + ^) (20b) 

In order to have non-zero potential minima above the surface of the top array, we have 
the following constraints 

c 2 B 0x > \B lx \ > 0, c c 2 B 0y > \B ly \ > 0, B 0x ^B 0y (21) 

The curvatures of the magnetic field at the centre of the traps and the trap 
frequencies (for the case of a harmonic potential) in the three directions are given by 

d 2 B d 2 B ld 2 B 47r 2 c 3| 



dx 2 dy 2 2 dz 2 a 2 



\B lx \ (22) 



^=. y = ^ = ^(^i^Y\B lx li (23) 
V2 a V m I 

The potential barrier heights in the three directions are given by 

AB X = AB y = c 4 \B lx \ (24a) 
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AB z =c 5 \B lx \ (246) 
C3, C4 and C5 are dimensionless constants which may be expressed in terms of Co 

C5 = (l+cg)* (25c) 
The barrier heights in the x-, y- and z-directions vary with z m i n as 

AS 21 = AB*c- fe ^"»"-*""™°) = AS y = ABjjfe - ***"**"-*™** ) (26a) 

AS 2 = AB z e- l(l """- z """°) (266) 

where 

A5£ = A_B^ = cia^e^ 2 """ (27a) 

AS 2 = Cl c 5 B 0x e- kz — (276) 

If Sis and Bi y are not subject to the constraint p6ajl we have AB X 7^ AB y . In some 
experiments, e.g., in studies of quantum tunnelling between lattice sites, it may be 
useful to be able to have different barrier heights along different axes. 

2.2.3. Three bias fields B\ x , B\ y and B\ z For the case of three bias fields B\ X) 
Bi y and B lz , the magnitude of the magnetic field is given by ifHfl) . The analytical 
expressions for B m i n , x m i n , y m in and z m i n in this case become very complex and we 
resort to numerical evaluation of to determine these quantities. Moreover, the 
condition for a symmetrical magnetic lattice in the case B z ^ imposes the constraints 

Box = Bq v , B\ x = Bi y (28) 

and under such conditions, B m i n = 0. Thus in the case of an infinite symmetrical 
magnetic lattice with three bias fields, B\ x , B\ y and Bi z , it is not possible to have 
non-zero potential minima. 

Nevertheless, we find that in the case of a finite magnetic lattice it is useful to be 
able to apply a bias field B\ z in order to compensate for the asymmetry introduced 
into the lattice by end-effects [21] associated with the finite number of magnets in the 
array (see section |3J). 



3. Numerical calculations for finite magnetic lattices 

To calculate magnetostatic potentials for arbitrary configurations, including finite 
periodic arrays of magnets of finite length and arbitrary cross section, the software 
package Radia [22] interfaced to Mathematica was used. Radia was also used to test 
various proposed configurations of permanent magnets that might support periodic 
arrays of microtraps with non-zero potential minima based on symmetry arguments. 
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3.1. Single finite periodic array of magnets with bias fields 

Figure (Hd) shows a numerical calculation of the magnetic field versus distance in the 
y-direction for the central region of a single finite array of n r = 1001 rectangular 
magnets, using the parameters a = 1 fim, t — 0.050 fim, l x — 1000.5 fmi and AttMz = 
3.8 kG (corresponding to the magnetization of our perpendicularly magnetized 
TbeGdioFesoCo^ magneto-optical films |13l and with no bias magnetic fields. 
Corrugations corresponding to the third-order spatial harmonic with period a/2 16 
[see (fTD|) and Ijlcll ] can be seen at distances very close (<C a/4ir) to the surface. At 
large distances, there are small residual corrugations with period a, which are due 
to 'end effects' [22 associated with the finite number (1001) of magnets in the array. 
This configuration may be used as a magnetic mirror for slowly moving atoms |14| . 

Figure^e) shows the effect of adding a bias field B\ y = —15 G. The magnitude 
of the magnetic field develops large corrugations with period a in the jy-direction, and 
2D magnetic traps with B m i n — appear in the potentials at z m i n — 0.540 /im. This 
configuration may be used as a spatial diffraction grating ^| 1181 1191 |2"U| for slowly 
moving atoms. 

Figure Elf) shows the effect of having bias fields B\ x — —20 G, B\ y = —15 G 
and B\ z = —0.09 G. The small value of B\ z was chosen in order to compensate for 
asymmetry introduced into the lattice by end-effects [2] associated with the finite 
number of magnets (1001) in the array. Non-zero potential minima with S m ,- n = 
\B\ X \ = 20 G appear in the potentials at z m i n = 0.540 fim. This configuration may 
be used as a ID periodic lattice of 2D magnetic traps or waveguides for slowly moving 
atoms. Relatively large values of B\ x and B\ y were used in this calculation so that the 
magnetic traps are sufficiently deep to be seen in the contour plot. In practice smaller 
values would normally be used. For example, if we use B lx = —2.73 G, B\ y = —0.22 G 
and Bi z = —0.09 G, non-zero potential minima with B m i n = \B\ X \ = 2.73 G appear 
in the potentials at z m i n — 1.216 /im. 

3.2. Two crossed layers of finite periodic arrays of magnets with bias fields 

Figure shows a numerical calculation of the magnetic field versus distance in the 
x-, y- and z-dircctions of a finite array consisting of two crossed layers of periodic 
arrays of parallel rectangular magnets. The parameters are n r = 1001, a — 1 /xm, 
s = 0.1 /im, t\ = 0.322 /xm (bottom layer), £2 = 0.083 [im (top layer), 4irM z = 3.8 kG, 
B lx = -4.08 G, B ly = -6.05 G and B lz = -0.69 G (see tabic [TJ. The values of 
Bi x , Bi y , s, t\ and ti satisfy the condition for a symmetrical infinite magnetic lattice 
with two bias fields [equations l|16a|l or (|166|l ]. while the value of B\ z was adjusted to 
compensate for asymmetry introduced into the lattice by 'end-effects' associated with 
the finite number (1001) of magnets in the array (see section 12. 2. 3(1 . The peaks at the 
ends of the array [figure|^a) and (b)] also arise from end-effects. The potential of the 
microtraps is close to harmonic in the region near the bottom of the traps [figure 
(g)-(i)]- 

Figure 0] shows a 3D plot and a contour plot of the magnetic field in the plane 
z = z m i n for the two crossed-layer structure of parallel magnets. The parameters used 
in this calculation (tabled column 3) lead to (n r — l) 2 = 10 6 magnetic microtraps with 
B m in — 2.73 G at the central minimum, which is located at d — 0.712 fj,m from the top 
surface, and magnetic field barrier heights AB X = B max — B min — AB V — 7.23 G and 
AB Z = 4.57 G. The potential minima in the outer regions of the array have slightly 
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Figure 3. Magnetic field produced by a periodic array consisting of two crossed 
layers of parallel rectangular magnets, (a) Near the edge in the x direction, (b) 
near the edge in the y-direction, and (c) in the z direction, for parameters given 
in table^ column 3. Curves (d)-(f) show the magnetic field in the central region 
of the lattice with the abscissa expanded. Curves (g)-(i) show the magnetic field 
in the central region of the lattice plotted (g) along a line (y = y m i„, z = z m i n ) 
parallel to the x-axis, (h) along a line (x = x m j n , z = z m i n ) parallel to the y-axis, 
and (i) along a line (x = x m in> V = Vmin) parallel to the z-axis. 




Figure 4. Magnetic field in the plane z = z m j n for a periodic array consisting of 
two crossed layers of parallel rectangular magnets, (a) 3D plot, and (b) contour 
plot, for parameters given in table Q column 3. 
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Table 2. Parameters calculated numerically using Radia (columns 3 and 6) and 
using the analytical formulae (columns 4 and 5) for the input parameters in tablelTl 
for (1) two crossed layers of periodic arrays of rectangular magnets (figure[3J, and 
(2) a single layer of square magnets with three thicknesses (figure 171. 

Config.(l) Config.(2) 

Analytical 



Parameter 


Definition 


Numerical 


B lz =0 


B u = -0.69 (G) 


Numerical 


&min (/^^) 


3? co-ordinate 
of potential 
minimum 


u.zou 


n 9^n 
u. zou 


n i 

U. lo4 


u. zou 


Vmin (tm) 


y co-ordinate 
of potential 
minimum 


0.250 


0.250 


0.294 


0.250 




z co-ordinate 
of potential 
minimum 


1.216 


1.216 


1.206 


0.952 


d (fim) 


Distance of 
central minimum 
from surface 


0.712 


0.712 


0.702 


0.732 




Magnetic field at 
potential minimum 


2.73 


2.73 


2.64 


1.10 


8 2 B / G \ 
dx^ cm 2 ' 


Curvature of B 
along x 


3.32 x 10 10 


3.31 x 10 10 


3.77 x 10 1() 


7.52 x 10 10 


d 2 B (_G\ 

dy^ * cm 2 ' 


Curvature of B 


3.32 x 10 10 


3.32 x 10 10 


3.77 x 10 10 


7.52 x 10 10 


along y 










ff (-St) 

as cm ' 


Curvature of B 


6.64 x 10 10 


6.64 x 10 10 


6.98 x 10 10 


1.50 x 10 11 




along 2 










AS 1 (G) 


Magnetic barrier 
height along x 


7.22 


7.22 


8.01 


8.10 


AB« (G) 


Magnetic barrier 
height along y 


7.23 


7.23 


7.78 


8.09 


AB 2 (G) 


Magnetic barrier 
height along z 


4.57 


4.57 


4.69 


5.45 



Table 3. Parameters calculated numerically using Radia for input parameters in 
tabic 111 for 87 Rb F = 2, mp = +2 atoms trapped in a magnetic lattice consisting 
of (1) two crossed layers of periodic arrays of rectangular magnets (figureEJ, and 
(2) a single layer of square magnets with three thicknesses (figure 171. 



Parameter 


Definition 


Configuration (1) 


Configuration (2) 


U min /k B (l*K) 


Energy of minimum of potential 


183 


74 


AU x /k B fjjiK) 


Potential barrier height along x 


485 


544 


Al/v/kB (pK) 


Potential barrier height along y 


486 


544 


AU*/k B (iiK) 


Potential barrier height along z 


307 


383 


u x /2n (kHz) 


Trap frequency along x 


232 


350 


uiy/2-K (kHz) 


Trap frequency along y 


233 


350 


<x) z /2tt (kHz) 


Trap frequency along z 


329 


494 


hw x /k B (jiK) 


Level spacing along x 


11 


17 


hujy/k B (flK) 


Level spacing along y 


11 


17 


HuJz/k B (jj.K) 


Level spacing along z 


16 


21 
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Figure 5. Distance of the central minimum from the surface d = z m i n — (ti + 
4 2 + s) (in units of a) versus bias fields B\ x , B\ y or B\ z , for the parameters given 
in table HI column 3. 




Figure 6. Potential barriers (a) and trap frequencies (b) in the x-, y- and in- 
directions, and B m i n (c) for the central minimum, as a function of d/a, for the 
parameters given in tabled column 3. 



different values for the barrier heights due to end-effects. The various quantities 
determined from this calculation are listed in table [5J column 3, along with values 
determined using the analytical expressions for an infinite symmetrical magnetic 
lattice with B\ z — (column 4) and B\ z — — 0.69 G (column 5). The values 
determined from the numerical calculations for the finite lattice with B\ z = —0.69 G 
are in excellent agreement with those determined from the analytical expressions for 
the infinite, symmetrical lattice with Bi z = 0. 

The height of the potential barrier in each direction i (i = x, y, z) for the trapped 
ultracold atoms is related to the magnetic field barrier height by 

A(7* = U; nax - U mm = m F 9 F ^AB' (29) 

The calculated parameters for 87 Rb atoms in the low magnetic field-seeking F = 2, 
rn F = +2 state are listed in column 3 of table The potential barrier heights 
are AU X = 485 /j,K, AU y = 486 \iK and AU Z = 307 fxK, and the trap frequencies 
are uj x = 2tt x 232 kHz, ui y = 2tt x 233 kHz and uj z — 2tt x 329 kHz. The trap 
frequencies may be scaled down, if necessary, by reducing cq = Bq x /Bq v or B\ x and 
Bi y [equation|J53J]. 

Figure 5 shows how the distance of the central minimum from the surface, 
d = z m i n — (t\ + t-x + s), expressed in units of the period a, decreases with increasing 
strength of the bias magnetic fields B\ x and B\ y . Figure Eta), (b) and (c) show the 
exponential increase in magnetic field minimum, potential barrier heights and trap 
frequencies with decreasing distance d of the minimum from the surface. 
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Figure 7. Periodic array consisting of one layer of square magnets with three 
different thicknesses. In (c) the locations of the central minima in the x-y plane 
are shown for a symmetrical magnetic lattice with parameters given in table Q] 
column 4. 



Figures 5 and El illustrate how the magnetic field minimum, potential barrier 
heights and the trap frequencies of the microtraps in the magnetic lattice may be 
controlled by varying the bias magnetic fields B\ x and B\ y to move the microtraps 
closer to or further from the magnetic array. We note that B\ x and B\ y should be 
varied simultaneously, according to l|16afl . while B lz is constant. 

3.3. Single layer of periodic arrays of permanent magnets 

A second configuration (figure 01 , which leads to qualitatively similar 2D magnetic 
lattices to the two crossed layers of parallel rectangular magnets, but which may be 
easier to fabricate, consists of a single layer of square-shaped magnets having three 
different thicknesses plus bias fields in the x-, y- and z-directions, where the thickness 
of the thickest magnet £3 = t\ +t2- The parameters used in the numerical calculation 
are listed in column 4 of table and the quantities determined from this calculation 
are summarized in the final columns of tables |21 and [3J 

4. Discussion and summary 

Using analytical expressions and numerical calculations we have shown that periodic 
arrays of permanent magnetic films plus bias magnetic fields can lead to 2D or ID 
magnetic lattices of microtraps having non-zero potential minima and controllable trap 
depth. Two configurations have been found that lead to 2D magnetic lattices with 
non-zero potential minima: the first consists of two crossed layers of periodic arrays of 
parallel rectangular magnets plus bias fields and the second consists of a single layer 
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of square-shaped magnets having three different thicknesses plus bias fields. These 
configurations lead to a symmetrical magnetic lattice with equal barrier heights in the 
x- and y-directions in the plane of the array if the bias fields B\ x and B\ y maintain 
a fixed relationship [given by (|16afl or (|16bt for the crossed array configuration] and 
if the bias field B\ z normal to the array is adjusted to compensate for asymmetry 
introduced by the finite size of the array. In some experiments, e.g., in studies of 
quantum tunnelling between lattice sites, it may be useful to be able to vary B\ x 
and B\ y independently in order to vary the relative barrier heights in the x- and 
y-directions. 

For arrays with micron-scale periodicity, the magnetic microtraps can have very 
large trap depths (~ 0.5 mK for the parameters chosen for the 2D lattice), allowing 
relatively warm atoms to be trapped, and very tight confinement. The barrier heights 
of the microtraps can be controlled by varying the bias fields B\ x and B\ v in the 
plane of the array, to move the traps either closer to or further from the surface. 
The numerical calculations for the crossed array configuration were performed for a 
1 mmx 1 mm magnetic lattice with period a = 1 /j,m, giving 10 6 lattice sites. It should 
be straight-forward to scale up the magnetic lattice, for example, to a 1 cm x 1 cm 
lattice with period a = 1 fim, giving 10 8 lattice sites. 

The permanent magnetic lattice configurations considered here should be suitable 
for trapping and manipulating small clouds of ultracold atoms prepared in low 
magnetic field-seeking states, including Bose-Einstein condensates and ultracold Fermi 
gases. A cloud of ultracold atoms could be loaded into the permanent magnetic lattice 
using, for example, a hybrid magnetic field structure comprising a current-carrying 
'U' quadrupole trap and a 'Z' Ioffe-Pritchard trap located beneath the permanent 
magnetic array on the atom chip JS] . Such a hybrid structure should allow ultracold 
atoms to be initially loaded from a mirror MOT into the 'U' surface MOT, then into 
the 'Z' magnetic trap, and finally into the 2D magnetic lattice. 

By using a small period (a ~ 1 \xm) and controlling the barrier height between 
the microtraps it may be possible to perform quantum tunnelling experiments such 
as the BEC supcrfluid to Mott insulator transition |3] in a 2D magnetic lattice. Some 
of the challenges will include the ability to fabricate permanent magnetic arrays 
with sufficiently smooth magnetic potentials and equivalent magnetic microtraps 
and to mimimisc the effects of the interaction of the ultracold atoms with the 
surface [23 EJ ESI EH1 EH EH EH] m order to preserve quantum coherence of the 
atoms in the magnetic lattice. For the magnetic lattices considered here, the potential 
minima are located about 1 fxm from the surface. At such distances the Casimir Polder 
force I2H] can be significant leading to an attractive component that lowers the 
barrier height, and losses due to thermally induced spin flips caused by interaction 
with the ambient temperature surface [23 EE1 E3 EH1 EHj 

can be important. However, 

spin-flip losses can be minimised by using magnetic films whose thickness (t < 0.4 /im 
for the magnetic lattices considered here) is much less than the skin depth and by 
use of suitable dielectric substrates with low electrical conductivity j2SJ. It should 
be possible to further reduce interactions with the surface, if necessary, by moving 
the potential minima further from the surface by decreasing the bias magnetic fields 
and/or marginally increasing the period of the magnetic lattice. Use of the BEC to 
Mott insulator transition for ultracold atoms trapped in 2D magnetic lattices could 
allow the preparation of a single qubit atom on each magnetic lattice site, which is 
important for scalable quantum information processing. 
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